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ABSTRACT 

We present the measurement of the two-point cross-correlation function (CCF) of 8, 198 Sloan Digital Sky 
Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Os- 
cillation Spectroscopic Survey (BOSS) at redshift z ~ 0.5 (0.3 < z < 0.9). The cross-correlation function 
can be reasonably well fit by a power-law model £,Qc(r) = (r/Vo)~ 7 on projected scales of r p = 2-25/r'Mpc 
with ro = 6.61 ±0.25/r'Mpc and 7 = 1.69±0.07. We estimate a quasar linear bias of £>g = 1.38 ±0.10 at 
(z) = 0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of 
~4x 1O i2 /i _i M , compared to ~ 10 13 /i -1 M Q characteristic host halo mass for CMASS galaxies. Based on 
the clustering measurements, most quasars at z ~ 0.5 are not the descendants of their higher luminosity coun- 
terparts at higher redshift, which would have evolved into more massive and more biased systems at low red- 
shift. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be 
dbg/dlogL = 0.20 ±0.34 or 0.1 1 ±0.32 (depending on different luminosity divisions) for quasar luminosities 
-23.5 > M,(z = 2) > -25.5, implying a weak luminosity dependence of quasar clustering for the bright end of 
the quasar population at z ~ 0.5. We compare our measurements with theoretical predictions, Halo Occupation 
Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of 
host halos, and the host halo mass distributions significantly overlap with each other for quasars at different 
luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find 
that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF 
equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These 
results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD 
approach and the need for additional information in populating quasars in dark matter halos with HOD. 
Keywords: black hole physics — cosmology: observations — galaxies: active — large-scale structure of Uni- 
verse — quasars: general — surveys 
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1. INTRODUCTION 

Quasars are powered by mass accretion onto supermas- 
sive black holes (SMBHs) at the center of massive galax- 
ies. Like galaxies, quasars are luminous tracers of the un- 
derlying dark matter, and can be used to map the large- 
scale structure of the Universe. Over the past decade, 
quasar clustering has been measured for large statistical sam- 
ples drawn from dedicated survey s, most no tably the Sloan 
Digital Sky Survey (SDSS, lYork et al l 120001) an d the 2dF 
QSO Redshift Survey (2QZ. ICroom et all 120041) . Build- 
ing o n earlier stud ies on small and heterogenous samples 
(e.g.. lShaverj|1984l) . the auto-correlation function of quasars 
has been measured with unprecedented precision for a wide 
redshift range (from z ~ 0.4 to z ~ 4) and a variety of 
quasar properties ("e.g.. [Porciani et al.l 120041 ICroomet al.l 
20051: iPorciani & Norbergll2006l iMyers et al.ll2006[ I2007albt 
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, 2008, 2009; da Angela et al 
l2009t llvashchenko et alJl2010HWhite etai]|2012h . 
and has been extended to the small-scale regime (5, 1 /r'Mpc , 
e.gJHennawi et al.l2 006: My ers et al.l2008tlShen et al.ll2010fc 
iKayo & Oguri 1 12012b . The clustering measurements have 
also been performed for Active Galactic N uclei (AGNs) se- 
lected at non-optical wavelengths (e.g.. IWake etalJ 1 2008 1: 
Gilli et alJ 120091: ICoil et alj 120091: iHickox et al l 120091 120TTI: 
Dono so et a J 120101: Cappellu ti et al.l 1201 Ot iKrumpe et alJ 
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l2Ol0l |20T1 iMivaii et al.ll2QTT1: lAllevato etafll20llh . These 
quasar/AGN clustering measurements revealed that quasars 
live in massive (~ 10 12 - 10 13 h~ l M Q ) dark matter halos, and 
constraints on the duty cycle of quasar activity can be in- 
ferred from the relative abundance of quasars and their host 
halos (e.g., ICole & Kaiserj|1989t iMartini & Weinberg! 120011; 
IHaiman & Huill2001l) . 

With quasar samples increasing in size, several attempts 
have been made to measure quasar clustering as a function 
of quasar luminosity. More massive halos are formed in rarer 
peaks of the de nsity fluctuation field and are more strongly 
clustered (e.g.. iBardeen et al.1 [19861 ICole & Kaiseri 11989b 
IMo & White] Il996t iSheth et alj I200lb . Gal axy clustering 



show s a strong dependence on luminosity (e.g., iN prberg et al. 



2001 


; Zehavi et al.l2005ll201 ltlCoil et al.l2006tlCoupon et all 


2012 


), indicating a good correlation between host halo mass 



and galaxy luminosity. On the other hand, quasar clustering 
studies to da te have failed to detect a strong luminosity depen- 
dence (e.g.. lAde lberger & Steidel 2005; Porciani & Norbergl 
2006 1 lMrors1st~ai1l2007ab Ida Angela et al.ll2008l; IWhite et al.l 
20121) . although She n et al.l (120091) reported a 2a detection for 
the most luminous quasars in SDSS Data Release 5 (DR5) at 
(z)~1.5. 

A weak dependence of quasar clustering on luminosity is 
expected if quasar luminosity is not tightly correlated with 
halo mass. Scatter between the instantaneous quasar lu- 
minosity and host halo mass dilutes any luminosity depen- 
dence of the clustering. Several semi-analytical cosmologi- 
cal quasar models have been constructed to make predictions 
broadly consistent with current constrai nts on the luminosity 
dependence of quasar clustering (e.g.jLidz et al.ll2006b IShenl 
2009; Shankar et aLllMOal: iConrov & Whitdl2013l" for recent 
work); more sophisticated approaches with dark matter-only 
simulations+semi-analytical galaxy form a tion models (e.g., 
Bonoli et alJ 12009b iFanidakis et"atl 120121: iHirschmann et al] 
20121). or with fully hydrodynamic cosmological simulations 
(g-g jThacker et al.ll2 009; Deg raf et al.lfeOl it IChatteriee et all 
120121) are underway. Precise measurements of the luminos- 
ity dependence of quasar clustering are important in quantify- 
ing th e scatter between quasar luminosity and host halo mass 
(e.g.. IWhite eTaTl 12008b IShankar et al.1 120 lOal) . which can in 
turn provide useful constraints on the correlation between 
black hole mass and halo mass, and on quasar light curve 
models (e.g lYu & Lul200l 120081; llfopkins et alj2005ll2008b 
[Sh^l2009[ICrotonll2009l;ICaoll2010l:IShanks et al.ll20 1 ltK 

The sparseness of quasars makes the measurements of the 
luminosity dependence of quasar clustering a nontrivial task. 
Fine bins in luminosity and redshift, while breaking the L-z 
degeneracy, lead to very noisy clustering measurements (e.g., 
Ida Angela et al.ll2008l) . hampering the d etection of a possible 
luminosity dependence. IShen et al.l (2009) used a flux-limited 
quasar sample covering a wide redshift range (0.4 < z < 2.5) 
in order to increase the statistics, but the resulting luminosity 
subsamples are mixtures over a range of quasar luminosity 
and redshift. 

One approach to mitigate such poor statistics is to cross- 
correlate the quasar sample with a much larger, galaxy sam- 
ple. On large scales, where linear bias applies, the cross- 
correlation function is determined by the auto-correlation 
functions of both sets of tracers. Using the cross-correlation 
technique, one can obtain a much better measurement of 
quasar clustering by boosting the pair counts, suppressing the 
shot noise from the small number of pairs in quasar auto- 



correlation measurements. In addition, the small-scale cross 
correlation between galaxies and quasars constrains the occu- 
pation of galaxies in quasar-hosting halos, and may hint on 
the triggering mechanism of quasar activity. 

There have been a number of studies on the cross correla- 
tion between galaxies and different types of quasars and Ac- 
tive Galactic Nuclei (AGN), i.e., optical-selected quasars, X- 
ray-, radio- and infrared-selected ( type 1 and typ e 2) AGNs 
(e.g. lAdelberger & SteideJ 12001 lUetall l2006t iCoil et al I 
120071 12009b IWake et al.l l2008t IP admanab han et alJ 120091: 
Donoso et al] 120 1 Ob iKrumpe et alJ 120101 120121 : iMiyaji et all 
1201 lb iHickox et al.1 12009T 1201 ll) . These studies generally 
found weak or no luminosity dependence of the large-scale 
quasar bias, although these measurements can be improved 
upon using larger samples. 

Here we use t he Tenth Dat a Rele a se (DR10), "CMASS" 
galaxy sample (IWhite etalJ 12011b lAnderson etaD 12012b 
Sanchez et al. 2012) fro m the Baryon Oscillation Spectro- 
scopic Survey (BO SS; [Schlegel, White & Eisen steinl 12009b 
iDawson et a"fll2012l) in SPSS I II jEisenstein et alj|201 ll) and 
the DR7 ([Abazaiian et al. 2009) sp ectro scopic quasar sample 
from SDSS I/II ( Schneider et al. 2010) to measure the cross 
correlation function of galaxies and quasars at 0.3 < z < 0.9 
((z) ~ 0.53). These samples represent the largest and most 
homogeneous spectroscopic samples to date for such cross 
correlation analyses, and enable us to derive one of the most 
stringent constraints on the luminosity dependence of large- 
scale quasar clustering in this redshift range. It also provides 
important clues on how galaxies and quasars occupy the same 
dark matter halos as functions of galaxy and quasar proper- 
ties, thus shedding light on the assembly process of quasars 
and their immediate environment. 

In this study we focus on the luminosity dependence of 
quasar linear bias at z ~ 0.5, although we also briefly touch 
on the occupation of quasars within dark matter halos. More 
detailed modeling and discussions on the other interesting as- 
pects of quasar-galaxy cross-correlation will be reported in 
future work. This paper is organized as follows: §|2]describes 
the quasar and galaxy samples used; the cross correlation 
function measurements are presented in §|3] we present a de- 
tailed discussion on our r esult s in terms of comparisons to the- 
oretical quasar mo dels ( 34.lt . Halo Occupation Distribution 
(HOD ) mo deling Q4.21 >. and mock catalog based interpreta- 
tion ( jj4.31 >: we conclude in $5] In the Appendix we present 
systematic checks of our correlation function measurements. 
Throughout the paper we adopt a flat ACDM cosmology with 
Q A = 0.726, h = 0.7 n h = 0.0457, er 8 = 0.8 and n s = 0.95 (e.g., 
Komatsu et alJ201 ll) . All errors quoted are 1 a statistical only, 
unless otherwise specified. Quasar luminosities are quoted in 
terms o f M,(z = 2), the absolu te i band magnitude normalized 
at z = 2 (iRichards et al.ll2006l) . 

2. THE DATA 

The SDSS I/II uses a dedicated 2.5-m wide-field tele- 
scope (iGunn et al.l |2006|) with a drift-sc an camera with 30 
2048 x 2048 CCDs dGunn et all 119981) to i mage the sky 
in five broad bands (ueriz: iFukugita et al.l 119961) . The 
imaging data are taken on d ark photometric nights of 
good se eing (|Hogg et aD 120011). are calibrat e d photomet- 
rically (ISmith et al l 12002b llvezic et al.1 12004b iTucker et al.1 
2006) and astrometricall y dPier et al.l 12003b. and object pa- 
rameters are_jnejisured (Luptonetal. 2001). Quasar can- 
didates dRichards et al.1 12002a!) for follow-up spectroscopy 
are selected from the imaging data using their colors, and 
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Table 1 

Summary of Quasar Subsamples. Nq G is the total number of quasar-galaxy pairs with r p < 50/r'Mpc and it < 70/r'Mpc in a given cross-correlation sample. 

The median redshift and magnitude are the pair-count (with r p < 50/r'Mpc and n < 70/r'Mpc) weighted median values of quasars. The last four columns list 
the best-fit power-law model correlation length of the CCF (with fixed slope 7 = 1.7), the galaxy linear bias, the linear bias of the cross-correlation sample fitted 
with the full covariance matrix and with diagonal elements of the covariance matrix. See g3]for details on subsamples and the estimation of correlation lengths 

and linear biases. 



# 


Sample 




Nq 


N G 


N QG 


2min 


Zmax 


^/.min 




<z> 


(Mi) 


D(7 = 1-7) 


b G 


£>QG 


, diag 
U 0G 





Full 






8198 


349608 


879352 


0.3000 


0.8999 


-28.693 


-22.576 


0.532 


-24.055 


6.614t-34 


2. 10 ±0.02 


1 70+ 006 


1.70 ±0.04 


1 


divl. 


_sl. 


_zl 


2726 


349608 


293098 


0.3003 


0.8998 


-25.115 


-22.576 


0.533 


-23.675 


6.682ttH 


2.11 ±0.02 


1 69+° 11 
1 °'-o. 1 1 


1.72 ±0.07 


2 


divl. 


_sl. 


_z2 


1075 


155888 


134524 


0.3003 


0.5320 


-23.819 


-22.576 


0.481 


-23.440 


6 390+ 0610 


2.03 ±0.04 


1 44+016 


1.42±0.10 


3 


divl. 


_sl. 


_z3 


1651 


193720 


135256 


0.5321 


0.8998 


-25.115 


-23.570 


0.589 


-23.942 


"""<••:::•.: 


2. 15 ±0.03 


1 90 +0 15 

1,u -0.16 


2.01 ±0.09 


4 


divl. 


_s2. 


_zl 


2738 


349608 


293640 


0.3002 


0.8999 


-25.541 


-22.808 


0.531 


-24.000 


6.841^1^ 


2. 10 ±0.02 


1 69 +008 


1.69 ±0.06 


5 


divl. 


_s2. 


_z2 


1068 


155888 


137808 


0.3002 


0.5319 


-24.171 


-22.808 


0.480 


-23.726 




2.03 ±0.04 


1 69 +01 ° 


1.68 ±0.08 


6 


divl. 


_s2. 


_z3 


1670 


193720 


133358 


0.5322 


0.8999 


-25.541 


-23.838 


0.591 


-24.294 


6.856^«i 


2. 15 ±0.03 


1 73+0-11 
'•' J -Q.12 


1.72 ±0.09 


7 


divl. 


_s3. 


_zl 


2734 


349608 


292614 


0.3000 


0.8993 


-28.693 


-23.208 


0.533 


-24.727 


f, 977+0.344 


2.11 ±0.02 




1.67 ±0.07 


8 


divl. 


_s3. 


_z2 


1069 


155888 


135812 


0.3000 


0.5319 


-26.851 


-23.208 


0.481 


-24.395 


6 823+ 571 

u -° -0.607 


2.03 ±0.04 




1.79 ±0.09 


9 


divl. 


_s3. 


_z3 


1665 


193720 


133933 


0.5327 


0.8993 


-28.693 


-24.204 


0.591 


-24.991 


5 303+ 533 


2. 15 ±0.03 


L58 -0.14 


1.52±0.10 


10 


divl. 


_s4. 


_zl 


837 


349608 


91081 


0.3004 


0.8993 


-28.693 


-23.915 


0.533 


-25.406 


6.804t«i 


2.11 ±0.02 


1 7Q+0.12 

'•''-0.13 


1.78±0.10 


1 1 


divl. 


_s4. 


_z2 


321 


155888 


41766 


0.3004 


0.5303 


-26.851 


-23.915 


0.482 


-25.043 


S.404^| 


2.03 ±0.04 


1 Q*J+0.18 
-0.20 


1.88±0.15 


12 


divl. 


_s4. 


_z3 


516 


193720 


42015 


0.5329 


0.8993 


-28.693 


-24.876 


0.592 


-25.622 


5.634*>|g 


2. 15 ±0.03 


1 ^Q+0.23 
'-■ Jy -0.2S 


1.42±0.17 


13 


div2. 


_sl. 


_zl 


2397 


249546 


283766 


0.3000 


0.5889 


-23.812 


-22.576 


0.484 


-23.564 




2.05 ±0.03 


1 67+°- 10 


1.63 ±0.07 


14 


div2. 


_sl. 


_z2 


1995 


78593 


136423 


0.3000 


0.4906 


-23.810 


-22.576 


0.448 


-23.420 


f, 7Q7+0.6I4 
°- ''-0.655 


2. 14 ±0.05 


1 55+0- 14 

'""-0.16 


1.52±0.10 


15 


div2. 


_sl. 


_z3 


402 


170953 


112867 


0.4907 


0.5889 


-23.812 


-23.369 


0.524 


-23.659 


6.429^' 


2.06 ±0.04 


1 69 +017 

l-tw_0 .19 


1.78±0.10 


16 


div2. 


_s2. 


_zl 


1443 


335123 


286117 


0.3005 


0.6980 


-24.315 


-23.812 


0.547 


-24.040 




2.11 ±0.02 


1 69+° 10 
1 °'-o. 1 1 


1.69 ±0.07 


17 


div2. 


_s2. 


_z2 


628 


178865 


123829 


0.3005 


0.5446 


-24.315 


-23.812 


0.499 


-24.018 


S 744+0-538 


2.05 ±0.03 


1 47+012 
'•^•'-0.13 


1.44 ±0.10 


18 


div2. 


_s2. 


_z3 


815 


156258 


132738 


0.5447 


0.6980 


-24.315 


-23.813 


0.592 


-24.066 


7 1 S0+°- 475 

'■ 1JU -0.499 


2. 12 ±0.03 


, 77+O.I6 
'■"-0.18 


1.84±0.10 


19 


div2. 


_s3. 


.^1 


4358 


349608 


306945 


0.3004 


0.8999 


-28.693 


-24.315 


0.578 


-24.741 


c q 71 +0.301 


2. 15 ±0.02 


1 7^+0.09 
-0.09 


1.74 ±0.06 


20 


div2. 


_s3. 


_z2 


624 


229499 


138601 


0.3004 


0.5747 


-26.851 


-24.315 


0.518 


-24.740 




2.09 ±0.03 


1 RS+013 


1.88±0.09 


21 


div2. 


_s3. 


_z3 


3734 


120109 


143922 


0.5748 


0.8999 


-28.693 


-24.316 


0.637 


-24.744 


6.259!^ 


2. 19 ±0.04 


1 77+O.O9 
'■"-0.09 


1.69 ±0.07 


22 


div2. 


_s4. 


_zl 


1966 


349608 


95949 


0.3019 


0.8999 


-28.693 


-25.000 


0.579 


-25.417 


6 030+ 05 ' 3 


2. 15 ±0.02 


1 75+0- 1 2 


1.70±0.11 


23 


div2. 


_s4. 


_z2 


188 


228104 


42244 


0.3019 


0.5738 


-26.851 


-25.003 


0.521 


-25.406 


5 936+ ™ 4 

J -"°-0.876 


2.09 ±0.03 


1 74+0.20 
'•'^-0.23 


1.69±0.17 


24 


div2. 


_s4. 


_z3 


1778 


121504 


45791 


0.5745 


0.8999 


-28.693 


-25.000 


0.644 


-25.419 


6 477+0.667 


2.20 ±0.04 


1 90+° 16 


1.84±0.13 



are arranged in spectroscopic plates dBlanton et al.l 120031) to 
be observed with a pair of fiber-fed double spectrographs 
dSmee et alj 12012b . The final (DR7) quasar ca talog from 
SDSS I/II was presented in lSchneider etal.1 d2010t) . 

T he BOSS survey is an ongoing program within SDSS 
III (lEisenstein et alJ 1201 lh . which is obtaining spectra for 
massive galaxy and quasar targets selected using photometry 
from SDSS I/II and new imaging data in the South Galac- 
tic Cap (SGC) in SDSS III. Targets are observed with an 
upgraded v ersion of the multi -object fiber spectrographs for 
SDSS I/II (ISmee et alJ l2012h . The BOSS spectra are re- 
duced and c l assifie d by an automatic pipeline described in 
iBolton et all (120121) . and the first p ublic data release of BOSS 
spectra is Data Release 9 (DR9) dAhn et all [20121) . In this 
work we use the unpublished Data Release 10 (DR10) for our 
galaxy sample, which contains BOSS spectra taken through 
July 2012, and surpasses the DR9 samples. 

2.1. Sample Construction 

We use the subset of qua sars in the SDSS DR7 quasar cat- 
alog (ISchneideretal.1120101). with UNIFORM _TARGET= 1 in 
the value-added catalog of IShen et alJ ([201 1). These quasars 
were uniformly targeted using the fin al quasar target selec- 
tion algorithm (Richards et al. 2002a) implemented in SDSS 
I/II, and cons titute a statistical sample suitable for clus tering 
studies (e.g.. IShen et al.ll2007l 120091: fRoss et al.l f2009). For 
the redshift range of interest here (z < 1), this quasar sample 
is flux limited to i= 19.1. The sky coverage of this uniform 
quasar sample is 6248 deg 2 . 

Two main galaxy samples are targeted in BOSS, with sep- 
arate color and magnitude cuts: the CMASS sample at (z) ~ 
0.55, and the LOWZ sample at z < 0.4. We choose the 
CMASS sample as our galaxy sample, as it has a larger red- 




Figure 1. Aitoff projection of the sky coverage of the cross-correlation sam- 
ples. The gray region shows the entire SDSS DR7 uniform quasar sample 
footprint, while the red region shows the current overlap with the DR10 
BOSS CMASS galaxy sample. 




Redshift 

Figure 2. Number density as a function of redshift for the DR7 uniform 
quasar and DR10 CMASS galaxy samples. We have limited both samples 
within 0.3 <z< 0.9. 
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Figure 3. Subsamples of quasars divided by quasar luminosity. The detailed sample definition is described in 32.2l and summarized in TablefT] The top panels 
show the distribution in the quasar luminosity-redshift plane, with different colors for the four different luminosity subsamples. Note that the red points overlap 
with the green points, i.e., the most luminous subsample is a subset of a less luminous subsample. The vertical dashed lines further split each luminosity subsample 
by the cross-pair-weighted median redshift. The bottom panels show the cross pair-weighted (with QG pair separations r p < 50A~'Mpc and n < 70/i~'Mpc) 
redshift distribution of quasars in each subsample (with the gray lines showing that for the full sample). The left and right columns are for Division 1 and Division 
2 in terms of quasar luminosity, respectively. 



shift overlap with our quasar sample. The total DR10 BOSS 
CMASS galaxy sample contains over 560 A: galaxies, which 
is approximately one half of the final BOSS CMASS galaxy 
sample. 

Since the CMASS galaxy sample has a narrow redshift dis- 
tribution that peaks around z ~ 0.55 and drops rapidly towards 
both ends, we have imposed a redshift cut, 0.3 < z < 0.9, to 
both the CMASS sample and the quasar sample. Fig.Q]shows 
the overlap between the CMASS galaxy sample and the DR7 
uniform quasar sample used in the current study, with a sky 
area of 4122deg 2 . Fig.|2]shows the redshift distributions of 
our final CMASS sample and quasar sample for subsequent 
cross-correlation analysis, with 349,608 galaxies and 8,198 
quasars in total. 

2.2. Quasar Luminosity Subsamples 

Since our primary goal is to investigate the luminosity de- 
pendence of quasar clustering, we divide our quasar sample 
into different subsamples by quasar luminosity. 

The redshift distributions of the quasars and CMASS galax- 



ies (e.g., Fig. |2} suggest that most of the pair contribution 
comes from a rather narrow redshift range around z ~ 0.5. 
Thus any redshift-dependent clustering is expected to be 
small. Nevertheless, we consider quasar subsamples divided 
by redshift-varying luminosity boundaries (Division 1), as 
well as by constant luminosity cuts (Division 2), as shown 
in Fig. [3] Division 1 enforces all subsamples to have the 
same redshift distribution, but the subsamples will overlap 
with each other in luminosity. Division 2 ensures there is 
no luminosity overlap in each subsample, but the effective 
redshift is slightly different for each subsample. We fur- 
ther split these luminosity subsamples by the pair-weighted 
quasar median redshift in each bin to create L — z subsam- 
ples, to investigate possible redshift evolution. Table Q] sum- 
marizes the luminosity and redshift boundaries and properties 
of these quasar subsamples. These redshift and luminosity 
boundaries were chosen to yield comparable pair counts for 
cross-correlation subsamples, except for the most luminous 
subsamples (divl_s4_* and div2_s4_*). 
We assign the effective luminosity and redshift to each 
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quasar subsample using the pair-weighted median values of 
quasar luminosity and redshift. 

2.3. Correcting for Fiber Collisions 

Due to restrictions of fiber placement during the BOSS sur- 
vey, two targets separated by less than 62" (corresponding 
to ~ 0.44/r'Mpc transverse comoving distance at z = 0.55) 
cannot be observed simultaneously on the same plate (tile), 
but can be both observed on overlapping plates. The BOSS 
tiling procedure uses optimized algorithms to maximize the 
number of galaxy targets in tile overlap regions, but there are 
still ~ 10% CMASS galaxy targets that do not have a spec- 
troscopic observation and are lost from the spectroscopic cat- 
alog. This fiber collision effect reduces the number of pairs 
on small (one-halo) scales and therefore lowers the clustering 
strength over these small scales. There are several schemes 
to compensate for the preferential loss of quasar-galaxy pairs 
due to fiber collisions: upweighting th e nearest spectroscopi c 
galaxies that have a collided target (lAnderson et al.l 120121) : 
assigning the photometric tar gets a redshift from the nearest 
spectroscopic neighbor (e.g., Zehavi et al. 2005); or using an 
algorithm that tracks the tiling g eometry and recov ers the true 
small-scale correlation strength dGuo et al.ll2012al) . 

Here we decided to use the upweighting scheme to re- 
cover the small-scale cross-correlation signal. In the case of 
our cross-correlation study, the spectroscopic observations of 
BOSS galaxies are completely independent of the spectro- 
scopic observations of the low-z SDSS-I/II quasars 16 , as the 
BOSS survey never plac es a fiber on a known low-redshift 
quasar dRoss et alJ 12012b . The upweighting scheme is thus 
equivalent to the nearest neighbor scheme such that both 
methods provide the maximum compensation for pair loss due 
to fiber collision. The information on the galaxy weights for 
fiber-collision (and a smaller fraction due to redshift failures) 
corrections is taken from the DR10 CMASS sample. 

2.4. Random Catalogs, Correlation Function Estimators, 
and Error Estimation 

We generate random catalogs for the CMASS galaxy sam- 
ple with the same angular geometry and redshift distribution 
as the data. The spectroscopic completeness f s (i.e., fraction 
of tar gets with fibers as signed) is a function of sectors (see 
e.g- lBlanton et a l. 2003, for the definition of sectors), and is 
taken into account by upweighting the galaxy points during 
pair counting. We already account for fiber collisions, so the 
spectroscopic completeness here does not include objects lost 
to fiber collisions. 

We estimate the ID and 2D redshift space correlation 
functions £ s (s) and £ s (r p ,Tr) using the simple estimator 
dDavis&Peebleslll983L DP): QG/QR- 1, where QG and QR 
are the normalized numbers of quasar-galaxy and quasar- 
random pairs in each scale bin, s is the pair separation in red- 
shift space, and r p (tt) is the transverse (radial) separation in 
redshift space. We shall comment further on this choice be- 
low. To reduce the effects of redsh ift distortions, we use t he 
projected correlation function (e.g., Davis & Peebles 1983) 

poo 

w p (r p ) = 2 dn £ s (r ; „7r) . (1) 
Jo 

In practice we integrate £ s (r p , 7r) to 7r max = 70/r'Mpc, where 

16 This situation is different from the cross-correlation between galaxies 
and quasars from the SDSS-I/II survey, where there is fiber collision between 
quasar targets and galaxy targets. 



Table 2 

Measurements of the cross-correlation function w p for the full sample and 
subsamples. The second column lists the total raw number of QG pairs in a 
given r p bin with it < 70/i~'Mpc, which can be used as a rough estimate of 

the robustness of the sample statistics. The last column lists the diagonal 
errors of the w p measurements, and the normalized covariance matrices are 

provided in Table|5] A portion is shown here for its content. The table is 
available in its entirety in the electronic version of this paper. 



sample 


r P 
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&w p ,diag 
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(ft-'Mpc) 




(ft-'Mpc) 


(ft-'Mpc) 





0.1155 


12 


2061.9628 


2440.1567 




0.1540 


23 


513.5358 


145.3023 




0.2054 


38 


464.1206 


127.6868 



the result is already converged for the scales considered in this 
paper. This upper-limit of 7T max will be taken into account in 
our subsequent modeling. For our fiducial £ s (rp,7r) grid we 
use a logarithmic binning in r p with Alogr ; , = 0.125 starting 
from r p m i n = 0.1 /r'Mpc and a linear binning in tt with Att = 
5r ! Mpc. 

There are different methods to estimate the statistical er- 
rors of the correlation function measurement, either inter- 
nally using bootstrap or jackknife resampling, or externally 
using mock catalogs (for a discussion, see, e.g jNorberg et alJ 
2009). Here we a dopt the jackknife resampling method (as 
was done in, e. g., iScranton et al.l 120021 iZehavi et"ail I2005t 
IShen et aT1l2007l) : we divide the clustering samples into jVj ac k 
spatially contiguous regions with equal area, and create Nj ac k 
jackknife samples by excluding each of these regions in 
turn. We create our jackknife samples using the pixelization 
sche me of STOMP 17 , wh ich has been used in other studies 
(e.g. JMcBridelt^l20Tlb . We measure the correlation func- 
tion for each of these jackknife samples, and the covariance 
error matrix is estimated as: 

Cov(/, j) = - 5x3 - ij) , (2) 

"jack ;=J 

where indices J and j run over all bins in the correlation 
function, and £ is the mean value of the statistic £ over the 
jackknife samples. The covariance matrix is generally domi- 
nated by the diagonal elements except for the large-scale bins, 
where correlations between adjacent £ bins become important 
due to common objects in these bins. 

We settled on 50 jackknife samples to estimate the covari- 
ance matrix. The normalized covariance matrix (also known 
as the correlation matrix) is defined as: 

CovO',;') 

Cov norm (i,7) = , (3) 

OiOj 

where of = Cov(/, z) is the diagonal element of the covariance 
matrix. By default we will use the full covariance matrix in 
our model fitting unless otherwise stated. Further discussions 
on error estimations and jackknife sampling are presented in 
the appendix. 

3. THE CROSS CORRELATION FUNCTION 

3.1. The whole quasar sample 

We show the projected correlation function w p for the full 
quasar and CMASS galaxy samples in Fig. [4] and tabulate 
the measurements in Table [2] Much of our focus will be on 

17 http : //code . google . com/p/astro-stomp/ 
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Figure 4. Projected cross-correlation function for the full quasar and 
CMASS galaxy cross-correlation sample. The black and cyan lines are the 
best-fit power-law model for the scale range r p =2 — 25 h Mpc with flexible 
power-law index 7 and fixed index 7fl x = 1.7. The red line is the best fit lin- 
ear bias model (i.e., the linear matter correlation function scaled by a constant 
bias) for the fitting range r p = 4— 16/r'Mpc. All fits were performed using 
the full covariance matrix. 



Table 3 

Quasar linear bias derived from bgQ and ba . The error bars are simply 
propagated from bqQ and be neglecting covariance. We only tabulated the 
results for the luminosity subsamples (e.g., the results for the L—z 
subsamples are too noisy to be useful). Note that the data for the most 
luminous subsample (s4) are a subset of the less luminous subsample (s3), 
so the bias measurements in these two bins are not independent. 
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0.547 
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0.578 


-24.741 


1.42 ±0.15 


div2. 


_s4. 


.zl 


0.579 


-25.417 


1.42 ±0.20 



the larger scales measurements, but it can be seen that we 
have a good detection of clustering to quite small scales. In 
particular, there are 842 QG pairs within r p < 1 /r'Mpc and 
7r < 70 /r'Mpc, allowing a fair estimate of the small-scale 
(one-halo) cross-correlation. 

We fit the measured CCF with a power-law model £(r) = 
(r/Vo)~ 7 over the projected scales 2 < r p < 25 /i~'Mpc to 
quantify the clustering strength on intermediate-scales. We 
can also estimate a linear bias bQc, i.e., 



(4) 



where w 



p, matter 



is the correlation function of the underlying 
matter at the redshift of interest, and bg C w bgbc where /?q 
and b G are the linear biases for the quasar and CMASS sam- 
ples respectively. 

To estimate the linear bias bQc, we use the linear matter 
correlat ion function com p uted using the linear power spec- 
trum in [Eisenstein & Hu ( 1999) under the adopted cosmol- 



ogy, estimated at the pair-weighted median redshift of the 
cross-correlation samples. Our investigations using mock cat- 
alogs (see ^4.31 l show that on scales r p < 4 /r'Mpc nonlinear 
and one-halo effects start to affect the linear bias, while at 
r p > 15 /r'Mpc residual redshift space distortion (RSD) ef- 
fects start to become important. Thus we narrow the fitting 
range to r p = [4, 16] /r'Mpc to estimate the linear bias, where 
a scale-independent linear bias seems to be a good approxi- 
mation (within 10%). Although we lose statistical power by 
excluding data points (i.e., only 5 bins of scale are used in the 
fitting), this procedure is preferred to avoid scales where non- 
linear effects, scale-dependent bias, and RSDs may affect the 
linear bias estimate. Nevertheless we tested varying r p bound- 
aries within [1 , 50] /i~'Mpc in the fitting and found all derived 
bQc values are consistent within 1 a, thus our estimate of bQc 
is robust against this detail. 

The correlation function is well fitted by a power-law model 
with ro = 6.61 ± 0.25 and 7=1 .69 ± 0.07 over the scales of 
2 < r p < 25 /r'Mpc (x 2 / dof = 6.54/7). On smaller scales, 
the correlation function significantly deviates from the best- 
fit power-law model derived from larger scales, and requires 
explicit modeling of the one-halo term. The fact that we detect 
significant clustering at r p < 1 /r'Mpc indicates that there are 
a population of satellite hosted quasars and CMASS galaxies 
in the cross-correlation sample (see discussions in §|4). 

The linear bias for the full cross-correlation sample from 
our simple fitting is /?qg = 1.70 ±0.06. In order to derive 
the quasar linear bias /?g we need to know the linear bias of 
CMASS galaxies be- For this purpose we have measured the 
auto correlation function (ACF) for the CMASS galaxy sam- 
ple using the standard DP estimator, and used the same fit- 
ting procedure to estimate be- However, we found that the 
best-fit be value does depend on the exact fitting range, given 
the substantially smaller statistical errors from the ACF mea- 
surement. To reduce the risk of contamination from small- 
scale non-linear clustering and large-scale redshift-space dis- 
tortion, we fit the CMASS ACF over the same scale range 
(r p = [4, 16] /r'Mpc) as for the CCF data, and derive b c = 
2. 10 ±0.02. Within this fitting range, the ratio of the CCF to 
the galaxy ACF is roughly constant, allowing use of the rela- 
tion bg C = bgbc to derive the quasar linear bias. The inferred 
quasar linear bias is /?g ~ 1.38 ± 0.10, consistent with the es- 
timated bQ ~ 1 .3 ±0.2 from the SPS S quasar auto-corr elation 
function measured at (z) ~ 0.5 (e.g.. lShen et alj |2069). This 
linear bias is also consistent w ith th e value derived using the 
HOD approach described in 84.21 and with the bias of the 
mock catalogs (which show a slight, slow decrease of the in- 
ferred bias from 4 /r'Mpc to 16/r'Mpc). 

Our derived CMASS galaxy bias value is somewhat larger 
than the estimated val ue of 1.8-2 in other ACF studies o f 
CMASS galaxies (e.g.. IWhite et"aT1l20il INuza et aljfe om. 
but is consistent with that derived in iGuo et al.l ( 2012H) 
based on the DR9 CMASS sample. This result is at least 
partly caused by the different methodology in estimating 
the bias. We also compared our ACF meas urement directly 
with those repo rted in other studi es (e.g ., IWhite et ail 120 lit 
lAnderson et ail I20H INuza et alj 120125 ; our measurement 
is systematically higher by ~ 10% over r p = 4- 16 /r'Mpc 
scales. To resolve this discrepancy we performed extensive 
tests upon our galaxy sample and the samples used in other 
studies, and found that this systematic difference is largely 
due to the usage of additional galaxy weights in the other 
studies. While there are good reasons to use those weights 
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Figure 5. Projected cross-correlation function for the quasar luminosity subsamples with the two luminosity divisions (see Fig. [3). The data points are mea- 
surements for that bin, with green symbols (within 2 < r p < 25/r'Mpc) indicating those used in the power-law model fitting. The w p data for the full sample 
is shown in dotted lines as a reference. The black dashed lines are the power-law fit to the fitting range r p = 2 — 25/r'Mpc with fixed slope 7 = 1.7, and the red 
lines are the linear matter correlation function scaled by the best-fit linear bias bqQ over the fitting range r p = 4- 16/r'Mpc. The sample number is marked in 
each panel (see Table[T]for sample information). 
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Figure 6. Projected cross-correlation function for the L— z subsamples with the two luminosity divisions. In each panel, the black shaded region is the \a range 
of the biased linear matter correlation function derived from the best-fit to the luminosity subsample indicated by the black number. The cyan and red points are 
the results for the two L—z subsamples of each luminosity subsample. For clarity we only show the data points over the 2 < r p < 25/r'Mpc range. 
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in these studies, it is not clear that they are applicable to our 
cross-correlation measurements. On the other hand, we tested 
the difference of usin g the simple DP estima tor and the more 
robust Landy-Szalay (lLandy & Szalayl ll993. LS) estimator, 
and found that the DP estimator over-estimates w p by only 
< 2% below r p = 10/i _1 Mpc and by ~ 10% at r p ~ 40/r'Mpc, 
which means the difference caused by using the simple DP es- 
timator is negligible. In general the statistical errors tabulated 
in Table[JJare significantly smaller than the systematic uncer- 
tainties in the galaxy bias estimation. Nevertheless, regarding 
the detection of the luminosity dependence of quasar bias, the 
exact value of the galaxy bias is not critical. 

3.2. Quasar subsamples divided in luminosity 

Fig. [5] shows the resulting cross-correlation function for 
each quasar luminosity subsample (i.e., no dividing in red- 
shift), and comparison with that for the full sample. For each 
luminosity subsample we show in Fig. [6] the results for the 
L— z subsamples. All the measurements are tabulated in Table 
2. 

Our current samples do not have a sufficient number of 
small-scale QG pairs (r p < l/r'Mpc) to probe the clustering 
difference on these one-halo scales when dividing our quasar 
sample in luminosity. To quantify the luminosity-dependence 
of the large-scale clustering strength we fit w p in the range 
of 2 < r p < 25/i _1 Mpc with the power-law model and in the 
range of 4 < r p < 16/r'Mpc with the linear-bias model. For 
the power-law model we fix the slope to be 7 = 1 .7, consistent 
with the best-fit slope for the full cross-correlation sample. 
The amplitude of the clustering is therefore measured by the 
best-fit correlation length ro and linear bias bQc- 

The best-fit values of ro and linear bias bQQ for different 
CCF subsamples are shown in Fig. [7] for the four quasar lu- 
minosity subsamples in each division. No significant dif- 
ference is detected among these subsamples. In Fig. [8] we 
present the w p values computed over wide linear r p bins with 
Ar p = 5 h Mpc for the four luminosity subsamples in the two 
divisions. These w p values represent the averaged correlation 
over these wide bins. Again, we see that while the value of w p 
depends on scale, there is no significant difference in cluster- 
ing strength between any of the samples on these scales. Our 
sample statistics are insufficient to probe potential luminosity 
dependence on r p < 1 /i _1 Mpc scales. 

One concern is that for Div 2 the effective redshift is slightly 
different for each luminosity subsample, and possible redshift 
evolution may complicate the interpretation. However, the 
difference in the linear growth factor over the probed redshift 
range (z ~ 0.45 -0.65) is only ~ 10%, and the evolution in the 
linear bias be of the CMASS galaxy sample over this redshift 
range is negligible (see Table [TJ. Thus the effect of redshift 
evolution is negligible for our samples, as expected, and we 
do not observe a significant difference when we further divide 
our luminosity subsamples in redshift (e.g., Fig.|9]l. 

4. DISCUSSION 

The improved measurement of quasar large-scale cluster- 
ing at z ~ 0.5, and the inferred luminosity dependence of 
quasar bias, can be used to study the evolution of the global 
quasar population and to test cosmological quasar models 
while the small-scale cross-correlation probes the immediate 
neighborhood of quasars and may hint at the triggering mech- 
anism of quasars. Since the statistics on the small-scale cross- 
correlation in the present study are still not sufficient for de- 



tailed studies (see Fig. [5]), much of our following discussion 
will focus on the large-scale quasar bias and its luminosity 
dependence, although we do attempt to model the small-scale 
clustering for the full cross-correlation sample. 

Quasars reside in dark matter halos, and the redshift evolu- 
tion of quasar bias can be used to understand the cosmic evo- 
lution of this population. A long-lived quasar population may 
passively evolve into their low er redshif t counterparts with a 
predi c ted bias evolution (e.g.. lFry||1996|; iTegmark & P eebles 
1998; Mo & White! [19961: iWhite et all l2007t Hopkins et al. 
2007a), and ca n be confronted with the observed quasar bias 
evolution (see 34. It . 

The observed luminosity dependence of quasar bias con- 
strains how well quasar luminosity correlates with halo mass. 
In a physical galaxy formation scheme, there are various cor- 
relations among halo, galaxy and BH properties such that a 
chain of L qso o M Bn o M gal o M h may form. If the BH 
mass is more directly connected to halo mass than to galaxy 
mass, we expect a simpler version, L qso «-» Mbh *h> Mh. In the 
simplest scenario, i.e., all quasars are shining at a constant Ed- 
dington ratio, and BH mass linearly correlates with halo mass 
with no scatter, we expect a strong luminosity dependence 
of quasar bias as a result of more luminous quasars living in 
more massive halos. In practice, there are inevitably curvature 
and scatter among these correlations, which will modify the 
resulting luminosity dependence of quasar bias. For instance, 
quasar luminosity at fixed BH mass may have a substantial 
dispersion, as a natural result from different fueling condi- 
tions; BH mass may not perfectly (and linearly) correlate with 
halo mass due to diversities in galaxy formation details. These 
scatters will produce a distribution of host halo mass at fixed 
quasar luminosity; the more these halo masses overlap in dif- 
ferent quasar luminosity bins, the less prominent will be the 
observed luminosity dependence of quasar bias. This effect 
will be further illustrated in the following discussion. 

4. 1 . Implications from large-scale clustering 

Fig. [10] presents the quasar/ AGN bias measured in differ- 
ent studies and comparisons to the bias of different galaxy 
samples. The three dotted lines show the bias of halos 
with c onstant halo mass Mh = 1,4,16 x 10 12 /z"'Mpc us- 
ing the iTinkere t al. (2005) halo bias formula 18 . The three 
dashed lines show the evolution of bias for a passive pop- 
ulation of tracers (e.g.. lFrv|[r9 96; Teg mark & Peebleslll998l 
IMo & Whitd [l996: Whit e et alJl2007HHopkins et alJl2007ah . 

These different samples probe different redshifts and lu- 
minosities, and are selected with different methods, thus a 
detailed comparison would be difficult. Furthermore, these 
studies used different methodologies to estimate the linear 
bias. Although in most cases the bias values derived with 
different methods agree to within 1 a, there are cases where 
they could differ sig nificantly (e.g. JPadm anabhan et al.l l2009t 
Krum peet alJl2012h . Keeping these caveats in mind, some 
general conclusions can be drawn from this figure: 

• Optically selected quasars appear to have a typ- 
ical halo mass b etween 10 12 - IP 13 h^Mp, (e.g. 
Croo m et al l 120051; iHopkins et al.l l2007at IShen et al.1 
20091: iShanks et alJl2'oT¥ over a wide redshift range. 
This result implies that most low-z quasars are not the 

18 Using alternative halo bias formula calibrated against simulations will 
yield sli ghtly differe nt results that are c onsistent within a factor of two (e.g., 
Sheth et al. 2001 ; Cohn & White 20®. 
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that the shape and amplitude of the cross-correlation function show no significant variation for different quasar luminosity subsamples. 



descendants of their high-z counterparts, which would 
have evolved into systems with relatively higher bias at 
low redshift. 

• There is no significant difference in the clustering 
strength between optical quasar samples and several X- 
ray selected AGN s amples at the same redshift (e.g., 
iKrumpe et alj 120121) . However, we note that these 
X-ray AGN samples only probe slightly fainter lu- 
minosities than the optical quasar samples, thus both 
types of active SMBHs are likely drawn from a simi- 
lar population, and therefore should trace a compara- 
ble halo mass range. There may be some hints that 
radio-selected AGNs have higher clustering than opti- 
cal quasars and X-ray selected AGNs (e.g..lWake et al.l 
l200llHickox et all2009tlDonoso et alJl2010h . 

• The galaxy populations from SDSS and BOSS are 
significantly more clustered than quasars/AGNs at the 



same redshift. By selection these galaxy samples are 
at the massive end of the galaxy population. Thus 
most low-z quasars are not shining within these massive 
galaxies. These massive galaxies may have experienced 
a brief quasar phase in the past to build up the central 
SMBH mass, and are therefore likely the descendants 
of high-z quasar host galaxies. 

At z ~ 0.5, the average stellar mass of the CMASS galaxy 
sample is ~ 2 x 10 n M Q (Maraston et al. 2012). This value 
corresponds to a black hole ma ss of ~ 4 x 10 & M Pl usin g the 
local M B h - Mbuige relation in iMarconi & Huntl (12003b and 
assuming all the stellar mass is in the bulge for CMASS 
galaxies. The average BH mass of the SDSS quasars is es- 
timated to be ~4x 10 7 Mq (assuming unity Eddington ra- 
tio) or ~ 3 x 1O 8 M (virial BH mass estimates from Shen et 
al. 2011). Since the SDSS quasars reside in halos that are 
typically a factor of a few less massive than CMASS galaxy 
hosts, either the quasar BH mass in these lower-mass galaxies 
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subsamples, connected by the dotted lines. No redshift difference is detected given the large error bars. 



is over-massive compared with the prediction from the local 
MBH-A^buiee relation, or the virial mass estimates for SDSS 
quasars are syst ematically overestimated (for the latter pos si- 
bility, see, e.g., She n et al.ll2008HShen & Kelly||20ial2"012l) . 

We now examine what constraints the luminosity- 
dependence of quasar bias at z ~ 0.5 can place on cosmolog- 
ical quasar models. First, we derive a quick constraint on the 
luminosity dependence of quasar bias by fitting a straight line 
to the data. For simplicity we neglect (small) correlated errors 
among these bias estimates due to the usage of the common 
galaxy sample in the cross-correlation measurements. Using 
the four luminosity subsamples in the two divisions, the slope 
constrained from the data is 

-^-=0.20 ±0.34 divl (5) 
dlogL 

= 0.11 ±0.32 div 2, (6) 

for -23.5 > Mj(z = 2) > -25.5. Thus the data are consistent 
with no luminosity dependence over this luminosity range. 

This weak luminosity dependence is in contrast to that 
of galaxy clustering (e.g., iNorberg et al.ll200Tl IZehavi et alj 
l2lxHl27mtrCoil et al.ll2008tlCoupon et alJl2012i) . The SDSS 
main galaxy sample at (z) ^0.1 shows a strong positive lumi - 
nosity dependence in galaxy clustering (IZehavi et al.11201 lb : 
b G (> L) x cts/0.8 = 1.06 + 0.21(L/L*) L12 , where L* corre- 



sponds to M r = -20.5. For the 0.4 < z < 0.6 galaxies in the 
Canada- France-Hawaii Tele scope Legacy Survey (CFHT-LS) 
sample (ICoupon et al.ll2012l) . b c (> L) = 1.166 + 0.288(L/L*) 
where L* corresponds to M* - 5 log/z = -19.81 (for all galax- 
ies). The luminosity dependence of galaxy bias for the CFHT- 
LS sample is shown in Fig. QT| and compared to that of the 
quasar bias derived in this work. We have assumed that the 
median quasar luminosity in our sample (M,(z = 2) = -24.055) 
corresponds to the galaxy threshold luminosity with the same 
bias, which incidently corresponds to a galaxy luminosity of 
w L* . Based on this comparison, a luminosity dependence of 
quasar clustering as strong as that for galaxies is ruled out at 
the ~ 95% (~ 2a) confidence level (CL). This result reflects a 
reasonably good correlation between galaxy luminosity (and 
stellar mass) and halo mass, a correlation that appears to be 
weaker between quasar luminosity and halo mass. 

The linear bias for a populati on of quasa rs at fixed luminos- 
ity L can be expressed as (e.g., Shen 2009b: 

f dP(M h \L) 
b Q (L)= / b h (Mh) \ j ' dM h , (7) 
J dM h 

where bh{Mh) is the linear bias of halos with mass M/,, and 
dP(Mh\L)j dMh is the distribution of host halo mass at fixed 
quasar luminosity L. If we define an effective halo mass 
(M h )(L) such that bh({M h )) = b Q (L), the dependence of (M h ) 
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passive population (e.g., Fry 1996; Mo & White 1996; Hopkins et al. 2008), started at three arbitrary high redshifts and matched to the measured linear bias of 
quasars at these redshifts. These biases derived in different work used different methods, and while th ey often agree within the reported error bars the re are cases 
when the reported error bars underestimate the systematic uncertainty in determining the bias (e.g., Padmanabhan et al. 2009; Krumpe et al. 2012), especially 
when the statistical uncertainty is small. With these caveats in mind, this figure suggests that quasars at different redshifts reside in halos with typical masses 
of a few lO 12 /t'Mq , and as such low-redshift quasars are not the descendants of their high-redshift counterparts, which would have evolved into more massive 
systems. The massive galaxies at z < 0.5 in the SDSS samples typically reside in ~ 10 13 /i~'Mq halos, and could be the descendants of z ~ 1 quasars. Right: 
Same as the left panel, but with the product of the linear bias and the linear growth factor D(z) as the y-axis. Thus constant large-scale clustering is denoted by 
horizontal lines in this plot. 

their BH mass correlates with halo mass as Mbh oc M^/ 3 
with no scatter (i.e., a "light bulb" model for quasars). The 
scalin g can be predicte d from some analytica l arguments 




(e.g., 



Figure 11. Comparison of the luminosity dependence of quasar bias de- 
rived in this work (symbols) with that of galaxies in the CFHT-LS sample 
(black solid line) at 0.4 <z < 0.6 I Coupon et al. 2012). We use open sym- 
bols for the second most luminous subsample (s3) in the two divisions to 
indicate the fact that it contains the most luminous subsample (s4). To map 
between quasar luminosity and galaxy luminosity we have assumed that the 
typical quasar luminosity in our sample (Mj(z = 2) = -24.055) corresponds 
to the galaxy luminosity with the same bias. Incidently we get a correspond- 
ing galaxy luminosity of ss L* . Note that the galaxy biases were derived 
for luminosity-threshold samples, and we have limited the galaxy luminosity 
within the range of 0.15-3L* , approximately the range probed by the CFHT- 
LS sample. The luminosity dependence of quasar bias is apparently weaker 
than that of the galaxy bias. 



on L determines the luminosity dependence of quasar bias. 
As a toy model, we parameterize a relation (M/,)(L) oc L a , 
A slope of a w 0.6 ~ 0.75 is consistent with a model in 
which all quasars are shining at fixed Eddington ratio, and 



Silk & Reesl [l998t IWvithe & Lo ebl 120031) or inferred 
from observation s of local dormant BHs (e.g.. iFerrare set2002t 
iBaes et aTll2003l) although scatter in the relation is expected. 
Any scatter in the M^-Mh relation, and dispersion in the Ed- 
dington ratio distribution, will lead to flattening in the (M/,)-L 
correlation (i.e., reducing a). Thus the level of observed lu- 
minosity dependence of quasar bias places a constraint on the 
scatter between halo mass and quasar luminosity for a given 
power-law slope in the intrinsic correlation. 

Fig. [12] (left) shows several realizations of this toy model 
with different values of a in dotted lines. Models with large 
a are less favorable compared with the data, although they 
cannot be completely ruled out given the uncertainties in the 
measurements. 

There are several more realistic, semi-analytical quasar 
models that can be confronted with this observational con- 
straint (see iUand Appendix B of White et al. 2012). It is be- 
yond the scope of this paper to compare these different mod- 
els in detail or use our measurements to constrain their model 
parameters (cf. Sha nkar et al.l l2010a b). 

As a simple demonstrati on, we c onsider one semi- 
analytical quasar model from lShenl (120091) . This cosmological 
quasar model assumes that quasars are triggered in halo ma- 
jor mergers, and adopts a quasar light curve model composed 
of an Eddington-limited accretion phase and a power-law de- 
caying phase. This model can reproduce a variety of quasar 
observables, including quasar clustering, luminosity function 
and Eddington ratio distributions over a wide redshift range. 
In Fig. [12] (left) we show the model predictions for the quasar 
bias as a function of luminosity at z = 0.5-0.6 as the gray 



12 



SHEN ET AL. 




-23.5 -24.0 -24.5 -25.0 -25.5 11.0 11.5 12.0 12.5 13.0 

M ; (z = 2) logM h (h _1 M ) 

Figure 12. Left: Comparisons between several model predictions and our measurement of the luminosity dependence of quasar large-scale linear bias. We use 
open symbols for the second most luminous subsample (s3) in the two divisions to indicate the fact that it c ontains the most luminous subsample (s4). For the 
dotted lines (i.e., power-law models with a = 0,0.3,0.6,0.75), the predictions are generated using the Tinker et al. (2005) halo bias for mula at z = .53, and 
normalized such that they are close to the measured bias for the full quasar sam ple. The gray band is t he prediction at z = 0.5 — 0.6 from the Shen (2009) model, 
and the blue dashed line is the prediction at z = 0.55 from the fiducial model in Conroy & White (2013, CW13) neglecting the satellite contribution (which serves 
to increase the bias in the fainter bins by about 5% while leaving the bright bins almost unchanged). Right: The distribution of host halo mass at fixed quasar 
luminosity from the Shen (2009) model, estimated at z = 0.5. 



shaded region. Although this model still predicts a mild in- 
crease in quasar bias with luminosity, it matches the data very 
well. The right panel of Fig. Q~2] displays the predicted distri- 
bution of halo mass for quasars at several fixed luminosities. 
There is considerable overlap in the range of halo masses for 
these quasar luminosities, which dilutes the bias difference of 
these quasars with different luminosities. The large dispersion 
in halo mass at fixed quasar luminosity is caused by both the 
scatter between halo mass and BH mass (or peak luminosity) 
and the lumin osity evolu tion of individual q uasars (see d iscus- 
sions in, e.g.Xidz et al.l l2006l: iWhite et al.ll2008l IShmll2009l: 
IShankaretal.ll2010al) . 

We also compare the data with the prediction from a simple 
model co nnecting halos and galax ies to quasars recently pro- 
posed bv lConrov & White! (12013b . This model is a "scattered 
light bulb" model which assumes a linear relation between 
galaxy mass and quasar BH mass, a lognormal distribution 
of quasar Eddington ratios, and a constant duty cycle. The 
free parameters in this model are tuned to match the observed 
quasar luminosity function over a wide redshift range. The 
predicted luminosity dependence of quasar bias at z = 0.55 
from their fiducial model (without satellite-hosted quasars) is 
shown as the blue dashed line in the left panel of Fig. [T2] 
This model predicts a luminosity de pendence tha t is slightly 
stronger than that predicted by the IShenl (|2009) model, al- 
though it is still consistent with the data within 1 a. Inclu- 
sion of satellite hosted quasars increases the predicted bias in 
the fainter bins by about 5% while negligibly changing the 
brighter bins. This marginally improves the agreement with 
our data. 

One might expect a stronger BH mass dependence of quasar 
clustering, because the additional scatter between the instan- 
taneous luminosity and BH mass (i.e., the Eddington ratio dis- 
tribution at fixed BH mass) has no effect here. Quasar BH 
masse s can be estimated with the v irial BH mass estimators 
(e.g., Vestergaard & Peterson 2006). We tested this hypoth- 
esis by dividing the quasar sample using virial BH masses 



estimated in IShen et all d201 ll) . but did not find any signifi- 
cant d ependence on virial BH mass (also see, e.g. JShen et alj 
2009). This result, however, could be due to the large statisti- 
cal and syst ematic uncertaint ies of these virial BH mass esti- 
mates (e.g.. lShen et al]l2008l) . or due to a large scatter in the 
intrinsic correlation between halo mass and quasar BH mass. 

4.2. Halo occupation distribution modeling 

Next, we attempt to model our CCF measurements with 
simple Halo Occupation Distribution (HO D) models (for a 
review on halo models, see, e.g. lCooray & Sheth 2002). This 
approach is an intuitive way to interpret the observed CCF, 
and can offer insights on how galaxies and quasars form in 
dark matter halos. 

We fix the galaxy HOD by ad opting parameters consistent 
with those in lWhite etafl d2011l) from modeling the CM ASS 
galaxy ACF, which reproduces our DR10 CMASS ACF mea- 
surement. The large-scale galaxy bias parameter from this set 
of HOD parameters is be = 2.00. For the quasar HOD, we 
focus on two types of parameterizations. Both types separate 
the contributions from central and satellite quasars 19 in halos, 
and they differ in the form of the central quasar HOD. In the 
first parameterization, the mean number of quasars located at 
the center of a halo of virial mass M is parameterized as 



(JVcenW) = 



1 



1+erf 



ZlogM-logM„ 



(8) 



This is a softened step function with characteristic mass scale 
M m j n and transition width of o i gM • We parameterize satellite 
quasars as a power law with a low mass rolloff, 



(AUM)) = exp 



Mo 
" M 



M 



(9) 



Such a quasar HOD parameterization is similar in form to 

19 In this work we use the term "satellite quasar" to refer to quasars hosted 
by satellite galaxies. 
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the galaxy HOD (e.g., IZheng etalJ 120051 120071) . and it is 
loosely motivated by cosmolog i cal hydrodynamic simu lation 
of AGN dDi Matteo et alJl200l iChatteriee et al.ll2012l) . This 
five-parameter model (M m ; n , tJ\ og M, Ma, M\, and a) has been 
applied to model the two-point auto -correlation functions o f 
(z) = 1 .4 and (z) =3.2 SDSS quasars dRichardson et al.l2012l) . 
The second quasar HOD parameterization adopts the same 
satellite HOD form, but it uses a log-normal form for the mean 
occupation function of central quasars, 



(N cen (M)) =/ cen exp 



(logM-logM cen ) 2 



(10) 



This parameterization has 6 parameters in total (3 for satellite 
HOD and 3 for central HOD). Compared to the 5-parameter 
model, it reduces the number of central quasars in massive ha- 
los. We will refer to the two types of HOD parameterizations 
as 5-par and 6-par models, respectively. 

For both parameterizations, we assume no correlation be- 
tween the occupation numbers of central and satellite quasars 
and between galaxies and quasars. We also assume that 
the spatial distributions of both quasars and galaxies in- 
side halos follow th e Navarro-Frenk- White (NFW) profile 
dNavarro et al.l U997). The variation and limitation of the 
quasar HOD parameterizations will be discussed after pre- 
senting the main modeling results. 

The calculation of the galaxy-quasar two-point C CF in 
the H O D framework foll ows sim ilar procedures in IZhengl 
d2004l) . IZehavi et all d2005h . and iTinker et al.l d2005l) . One 
improvement we have in the model is to incorporate the ef- 
fect of residual redshift-space distortion (RSD) when com- 
puting the project ed CCF from th e real-space CCF, by apply- 
ing the method of lKaiseri (119871) to decompose the CCF into 
monopole, quadrupqle, an d hexadecapole moments (also see 
Ivan den Bosch et aill2012t J. Tinker, private communication, 
2009), which improves the modeling on large scales as we 
will see later. 

We model the cross-correlation between CMASS galaxies 
and the full sample of quasars at the pair-weighted redshift z = 
0.53. We include the quasar number density in calculating x 2 , 
adopting a value of 2 x 10~ 6 /i 3 Mpc~ 3 with a 20% fractional 
error (see Figure 2). A Markov Chain Monte Carlo method is 
applied to probe the parameter space. 

The main results from the HOD modeling are shown in Fig- 
ureQj] In FigureQjja), the solid curve is the best-fit w p from 

the 5-par model, with % 2 /dof=26.6/18. The value of x 2 is 
about 1.4er higher than the expected mean value 18, which is 
mostly contributed by the three points between 20/i _I Mpc and 
40/r'Mpc. While it is an acceptable fit, the slightly higher 
X 2 may indicate that the model needs further improvement or 
that the error bars and covariances on large scales are under- 
estimated. The dashed curve shows the predicted w p with the 
best-fit HOD if the residual RSD is not included in the model. 
As expected, on scales much less than 7r max = 70/r'Mpc, the 
effect of residual RSD is small. However, on scales close to 
T^max, the effect starts to appear, e.g., about 40% lower in w p 
at r p ~ 50/r'Mpc if the residual RSD is neglected. The x 2 
from the w p with no RSD becomes x 2 /dof=33.3/18, clearly 
demonstrating that including the residual RSD does improve 
the fitting significantly. 

The best-fit mean occupation functions for the 5-par model 
are shown in Figure [T3l fe). which can also be interpreted as the 
mass-dependent duty cycle of the quasars in the full sample, 



i.e., the fraction of halos hosting active quasars in the full sam- 
ple. For central quasars, a large transition width of the soft- 
ened step function makes (N c&n (M)) behave like a power law 
with an index of ~ 0.8 above lO'V^'M©. Satellite quasars 
(with power law index ~ 1.07 in (N s - dt ) at the high mass end) 
start to dominate around 10 14 /! _I M©. The overall occupa- 
tion function resembles a power law with index ~ 0.95. The 
shaded regions delineate the envelopes from the first 68.3% of 
the models after sorting them in ascending order of x 2 , which 
give us some idea of the constraining power of the CCF on 
the quasar HOD. For central quasars, the high-mass end is not 
well constrained - the fast drop in halo mass function toward 
the massive end makes quasars in massive halos contribute lit- 
tle to the large scale bias and number density of quasars. For 
satellite quasars, the constraints are tighter around the mass 
scale where they become comparable in occupation number 
to the central quasars. This mass scale also corresponds to 
the mass range of halos that have a significant contribution to 
small-scale galaxy-quasar pairs. Other than this mass range, 
the constraints on satellite HOD are loose. 

Multiplying the best-fit mean occupation function with the 
differential halo mass function, we obtain the contribution to 
the quasar number density from halos of different masses, 
as shown in Figure IT3l e). With appropriate normalization, 
the curve also gives the probability distribution of the host 
halo mass of the quasars in the full sample. While peaked 
around 1O 12 /i _1 M , the host halos have a wide distribution 
in mass, about 4 dex in a full-width-half-maximum sense. 
Marginalized over all models, the median host halo masses 
for central and satellite quasars are logM me d, ce n = 1 1 -60^q 
and logM mec j. S at = 13.74+j} 2 g, respectively. 

Figure \l3\ e) demonstrates that satellite quasars (dashed 
curve) clearly make a non-negligible contribution to the full 
sample. The strong small-scale clustering in the data requires 
the existence of satellite quasars. Otherwise, the small-scale 
w p would become shallower. The satellite fraction marginal- 
ized over all models is / sa t = 0.068^] ( me tnm curve m Fig- 
ureptt/)). 

With the adopted HOD parameterization, the 5-par model 
successfully reproduces the observed galaxy-quasar CCF. The 
central quasar occupation function appears to be a signifi- 
cantly softened step function (eri og M = 2.73^y2i)- Such a large 
transition width implies a large scatter in quasar luminosity 
at any given halo mass. The large transition width also leads 
to a wide mass range of host halos, which even extends to a 
few times 10 9 /z _1 Mq, a regime for dwarf galaxies. This result 
of low mass halos does not appear to be reasonable. Could 
it be an artifact of the parameterization of the 5-par model? 
The (N cen (M)) function is parameterized to be monotonically 
increasing with mass towards an asymptotic value of unity 
(although it never reaches unity in the mass range of interest). 
There are only two free parameters in (N cen {M)}, making a 
relatively tight connection between the high-mass end and the 
low-mass end HOD. For example, while a higher (N cen (M)) 
at the high mass end helps to reproduce the small-scale clus- 
tering, it increases the large-scale bias, and as a response, the 
occupation function must extend to low-mass halos to reduce 
the large scale bias. 

The 6-par model can explore the parameterization limita- 
tion, which allows the high-mass occupation function of cen- 
tral quasars to cutoff exponentially. It tends to mimic the lack 
of quasar activity in high mass halos where gas accretion is 
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Figure 13. Results from HOD modeling of the cross-correlation between galaxies and the full sample of quasars. Panel (a): HOD fit to the projected galaxy- 
quasar CCF. The solid curve is the best-fit from the 5-par HOD model with the effect of residual redshift space distortion (RSD) included. The shaded region is 
the envelope of the fits from the 68.3% of the models with the smallest x 2 values in the MCMC chain. The dashed curve is the predicted w p with the above best-fit 
HOD, if the effects of residual RSD were not included. Panel (b): The best-fit mean occupation function of quasars (solid) from the 5-par model, decomposed 
into its central (dotted) and satellite (dashed) components. The red and blue shaded regions are envelopes from the 68.3% of models with the lowest \ 2 values 
for the central and satellite mean occupation functions. Panel (c): Same as (b), but from the 6-par model. Panel (d): The fraction of satellite quasars in the full 
sample derived from the HOD modeling. The thin and thick curves are from the 5-par and 6-par models, respectively. Dotted lines enclose the central 68.3% of 
each distribution. Panel (e): The contribution to the quasar number density as a function of halo mass, decomposed into central (dotted) and satellite (dashed) 
quasars, from the best-fit 5-par model. The curves are obtained from the product of the mean occupation functions and the differential halo mass function. The 
curves are also proportional to the probability distribution of host halo mass of quasars. Panel (f); Same as ( e), but from the 6-par model. See the text for details 
on the 5-par and 6-par models. 



likely suppressed. With this 6-par model, we find an almost 
equally good fit to w p , with x 2 /dof=26.1/17, and the best-fit 
curve is similar to that in Figure fT3l a). The constraints on the 
mean occupation functions (indicated by the shaded regions in 
FigureQjfc)) become less tight, especially for central quasars. 
The host halo mass for central quasars now has a much nar- 
rower distribution (see Figure [131 0). which is in a better 
agreement with the prediction from the Shen (2009) model 
(See the right panel in Fig.[T2l. Marginalized over all models, 
the median host halo masses for central and satellite quasars 
are logM med , cen = 11.85^;| and logM med , sat = 13.66^, re- 
spectively. The satellite fraction from the 6-par model is 
/ 8at = (see the thick curve in Figure El^)). 

The high satellite fraction from either model is a some- 
what surprising result. W ith a similar 5-par parameterization, 
iRichar dson et aQ (12012b model the 2-point auto-correlation 
function of 0.5 < z < 2.5 (z = 1.4) SDSS quasars and infer a 
satellite fraction of (7. 4 ± 1.3) x 10~ 4 . Also from HOD mod- 
eling of quasar clustering, Kavo & Oguri] (12012b infer a satel- 
lite fraction of 0.054+ogjg for 0.6 < z < 2.2 quasars. Although 
our result is close to the lat t er one , the parameterizations are 
different — Kayo & Oguri] (120121) assumes that both the cen- 
tral and satellite quasar occupation functions have the same 
Gaussian form, differing only in the amplitudes. The satel- 



lite fraction is mainly determined by the small-scale cluster- 
ing. In detail, for our quasar-galaxy CCF modeling, the result 
would depend on the assumptions about the correlation be- 
tween galaxies and quasars inside halos and about the spatial 
distribution of satellite quasars and galaxies inside halos. This 
again highlights the ambiguity in HOD parameterizations for 
the quasar population. 

One important distinction is that the quasar satellite frac- 
tion in our HOD model is not the fraction of binary quasars 
(quasar pairs on 1-halo scales). Many of the massive halos 
will only have one satellite quasar and no central quasar, thus 
the actual binary quasar fraction would be substantially lower 
than the satellite fraction. We still designate these quasars as 
satellite quasars (even though they are the only quasar in the 
halo) because they have a distinct intra-halo spatial distribu- 
tion compared to central quasars in our HOD modeling. 

The clustering measurement can be well fit using different 
HOD parameterization, as demonstrated by our 5-par and 6- 
par models. That is, there exist large degeneracies in quasar 
HOD from the clustering data alone. In addition to the 2- 
point correlation functions, we need other observables (e.g., 
pairwise velocity distribution) to break the degeneracies and 
constrain the connection between quasars and halos. We also 
need to rely on theoretical work for a more physically moti- 
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Figure 14. Top: the mean (total) occupation number of q uasar s and galax- 
ies for the two quasar HOD parameterization described in 34.21 The galaxy 
HOD is the CMASS HOD shifted to lower mass scales to mim i c a L > L* 
galaxy sample, which seems consistent with that in Coupon et al. (2012), and 
roughly matches the large-scale clustering of quasars. Bottom: the ratio be- 
tween the mean occupation numbers of quasars and galaxies. The shaded 
region indicates the 68.3% confidence range. For both quasar HOD parame- 
terizations the ratio of quasars to galaxies rises to a plateau at the high-mass 
end, but the uncertainties are too large to confirm or rule out a decline in the 
quasar fraction (per galaxy) in > 10 i4 Mq halos (e.g., clusters of galaxies). 



vated HOD parameterization to model quasar clustering. 

We also tried to model the HOD for our quasar luminosity 
subsamples, but the constraints are poor given the increas- 
ingly larger measurement uncertainties. Therefore we defer a 
more detailed HOD modeling of the luminosity dependence 
of quasar clustering to future work with improved clustering 
measurements (especially on small scales, see discussions in 
§ 4.4). The large-scale quasar bias for the full sample from 
our HOD modeling is: b = 1.27!°$ (5-par) and b = 1.26![J$ 
(6-par), which are slightly lower, but consistent with our esti- 
mation in ^3] within la. 

Finally, we comment on whether quasars are under- 
represented in massive halos by examining the ratio of quasars 
to galaxies as a function of halo mass. Fig. [T4l shows the ratio 
of (central+satellite) quasars to galaxies as a function of host 
halo mass, for the two HOD parameterizations above. For 
the galaxy HOD we have simply shifted the CMASS HOD 
to lower mass scales to approximate a L > L* g alaxy sample, 
which seems to be consistent with the results in Coupo net al.l 
d2012l) . and roughly matches the large-scale clustering of 
quasars (see Fig. [TT| and caption thereof). The quasar-to- 
galaxy ratio rises to a plateau at high halo masses in both 
HODs, but the uncertainties are large and we cannot con- 
firm or exclude a decline of quasar fraction (per galaxy) in 
> 10 14 M© halos (e.g., clusters of galaxies). 

We tabulate the best-fit quasar HOD parameters and the 
adopted CMASS galaxy HOD parameters in Table [4] but we 
caution that the quasar HODs are merely for future reference 
purposes and not for detailed physical interpretation, given 
the large degeneracies discussed above. 

4.3. Mock catalog based interpretation 

We now consider a mock catalo g based approach to in- 
terpret the observed CCF (e.g. Padmanabhan et al. 2009; 



Table 4 

The adopted CMASS galaxy HOD parameters and the best- fit parameters 
for the two quasar HOD parameterizations described in 34.21 All masses are 
in units of /t'Mq . We caution that the quasar HODs are merely for future 
reference purposes and not for detailed physical in terpr etation, given the 
large degeneracies discussed in 34.21 



CMASS HOD 
Eqs. (8) and (9) 



5-par quasar HOD 
Eqs. (8) and (9) 



6-par quasar HOD 
Eqs. (9) and (10) 
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I White et all 1201 U IConrov & White! [20 131) . Compared with 
analytic implementation of the HOD ( jj4.2l) . the mock-based 
approach directly uses simulated halo catalogs, thus avoid- 
ing using any specific fitting formulae for the halo bias and 
abundance. Unfortunately it can be subject to finite volume 
and finite resolution limitations. The basis of our catalogs is 
a 2048 3 particle N-body simulation of the ACDM cosmology 
in a 700/i~ 1 M pc box run with the TreePM code described in 
I White! d2002l) . This simulation has sufficient volume to probe 
the CCF on the scales of relevance here while retaining suf- 
ficient force and mass resolution to resolve the halos hosting 
CMASS galaxies and quasars. 

We can populate the halos in the simulation using differ- 
ent models for the relevant objects. The CMASS galaxies 
are plac ed in the halos using a HOD similar to that described 
in jj4.2l The paramet ers are adju s ted to fit the small-scale 
clustering measured in White et al. (201 1) and the large-scale 
clustering measured in And erson et al.l (120121) for CMASS 
galaxies. Since our purposes are primarily illustrative, we 
simply chose one model which provides a good fit without 
attempting to propagate the uncertainty in this model. This 
best-fit model is a very good fit to the data. For the quasars 
we chose two different models based on the framework in 
IConrov & Whltel J2013l CW13 for short). The CW13 frame- 
work assumes there is a linear relation between galaxy stellar 
mass and BH mass with a scatter, and that the BH shines as a 
quasar with a constant duty cycle, with its luminosity drawn 
from a lognormal distribution with a constant mean Eddington 
ratio. This simple model can reproduce the quasar luminos- 
ity function and large-scale quasar bias for a wide range of 
redshifts. 

For both quasar models we consider the cross-correlation 
on both large- and small-scales is independent of the overall 
duty cycle of the quasars — a random dilution of the sam- 
ple returns the same clustering on average. The first model 
assumes quasars live at the centers of dark matter halos with 
the quasar luminosity set by the stellar m ass of the galaxy 
most likely to be hosted by such a halo (as in Conroy & White! 
120131) . In the second model, quasars live in both central and 
satellite galaxies, with th e quasar luminosity s et by the stellar 
mass of the galaxy (as in lConrov & Whitel2013[) . Comparison 
between the two models shows the impact of quasars populat- 
ing satellite galaxies. 

Fig. Q3] shows the CCF comparisons of our mock predic- 
tions with the data, for the three luminosity subsamples: 13, 
16 and 19 in Division 2 (see Table Q]). In each panel, the 
black line with error bars is the measured CCF, and the red 
(CW13-cen) and cyan (CW13-all) points are our mock pre- 
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Figure 15. Compaiisons between the measured CCF and predictions from our mock catalogs, for the three luminosity subsamples 13, 16 and 19 (see Table[TJ. 
In each panel the black line with error bars is the measured CCF, the red open squares are the prediction for mock quasar model (1) and the cyan filled circles are 
the prediction for mock quasar model (2). The errors on the predicted CCF are smaller than the observational errors, and are suppressed for clarity. See text for 
details on the mock catalogs and interpretations. 

samples shown in Fig. [15] This satellite fraction is si milar to 
that inferred from the 6-par HOD model discussed in § 14.21 In 
reality, the situation may be more complicated such that cen- 
tral galaxies might be less likely to host a quasar than satellite 
galaxies in the most massive halos (e.g., clusters), which will 
lead to changes in the satellite fraction. In addition, just as 
for our HOD modeling, any enhanced probability of finding 
close galaxy-quasar pairs (e.g., if quasars are triggered during 
interactions with companion galaxies) will change our mock 
interpretation (which assumes galaxies and quasars are statis- 
tically independent when populating the halos). Additional 
observations of quasars in groups and clusters are required to 
probe these possibilities. 

For our mocks, the mean quasar occupation number and the 
distribution of host halo mass differ in detail from our best- 
fit HOD models in 34.21 which again highlights the fact that 
there is a broad range of HOD parameter space that can ac- 
commodate the observed CCF. 

A side product of our mock-based modeling is a prediction 
for the scale-dependence of the bias for the CCF. In Fig.[T6lwe 
show the ratios of the CCF of our mock catalogs to the auto- 
correlation function of the underlying dark matter computed 
from the linear and non-linear matter power spectra from our 
simulation. The linear bias is approximately constant over 
scales ~ 4- 16/r'Mpc. It is on the basis of this modeling that 
we have chosen the fitting range quoted in §[3] 
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Figure 16. Linear and non-linear biases of the CCF from one of our mock 
catalogs. The underlying matter correlation function was computed us- 
ing the linear and non-linear power spectra from the simulation directly. 
The shaded region encloses the ±5% range of the median non-linear bias 
within r p = 4- 16/r'Mpc. Both the linear and non-linear biases show scale- 
dependence. The non-linear bias is computed using the projected correlation 
function including redshift space distortions while the linear bias calculation 
does not include redshift space distortions. For scales 4 < r p < 16/r'Mpc, 
the linear bias is roughly scale-independent. This result motivated our choice 
of the fitting range in deriving the linear bias in i}5] for which the effects of 
scale-dependent bias and redshift space distortions are negligible. 

dictions for quasar model (1) and (2), respectively. Model (1) 
where quasars only populate central galaxies does not pro- 
vide a good match to the small-scale CCF. On the other hand, 
Model (2) where quasars populate both central and satellite 
galaxies provides a good match to the overall CCF for three 
luminosity subsamples (although the model may over-predict 
the CCF a little on scales of a few /r'Mpc for sample 19). 
The reason that the predicted CCF does not vary much over 
the three quasar luminosity bins is that there is substantial 
overlap in the host halo mass range for quasars in the three 
bins, due to the significant scatter between host galaxy stel- 
lar mass and instantaneous quasar luminosity in the CW13 
model (~ 0.4dex). Since in Model (2), quasars are randomly 
subsampled from galaxies regardless of their positions (with 
scatter), the overall satellite fraction of quasars is roughly the 
same as for galaxies, i.e., / sat ~ 10% for the three luminosity 



4.4. The future 

Given the weak luminosity dependence of quasar cluster- 
ing, one must considerably improve the errors on the mea- 
surements to firm up a detection. In addition, it is desirable to 
have a larger lever arm in quasar luminosity, since the change 
in quasar linear bias with luminosity is slow. With the cross- 
correlation technique the galaxy sample limits us to a fixed 
area of sky. To go brighter we need to work at the highest red- 
shift available (both because of volume effects and because of 
the z-dependence of the luminosity function). To go fainter 
we need to probe to dimmer objects in the same area of sky. 

A major discriminant between quasars models lies in the 
less luminous quasars (below L*). In older, or more simpli- 
fied, models these quasars arise from low-mass black holes 
accreting at close to the Eddington rate, whereas in most mod- 
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ern models a significant fraction of them arise from higher 
mass black holes accreting at a lower rate (and the prevalence 
of low accretion rate black holes is particularly pronounced in 
the redshift range of interest here). 

Unlike most galaxy clustering measurements (especially 
those from SDSS), quasar clustering measurements are still 
limited by statistical errors. Our current cross-correlation 
sample only includes ~ 2/3 of the final CMASS galaxy-DR7 
quasar overlap sample. Thus we expect some improvement 
in the clustering measurements using the final data release 
of BOSS. The signal-to-noise ratio for Poisson noise domi- 
nated regimes (e.g., at small scales) will increase by a factor of 
~ 1 .2. For large-scale bins where errors are correlated, we ex- 
pect improvements somewhat smaller than this. In any case, 
the final cross-correlation sample will have a more uniform 
sky coverage than the current sample, which may eliminate 
some systematic problems. 

5. CONCLUSIONS 

In this paper we presented the cross-correlation function 
measurements between quasars and galaxies at z ~ 0.5 using 
a spectroscopic quasar sample from SDSS DR7 and a BOSS 
CMASS galaxy sample from SDSS-III DR10. Our cross- 
correlation sample contains 8,198 quasars and 349,608 BOSS 
(CMASS) galaxies. Our main results are the following: 

• The CCF can be well described by a power-law model 
£,qg = ( r / r o)~ 7 for scales r p = [2,25]/r'Mpc with ro = 
6.61 ±0.25 A _1 Mpc and 7 = 1.69 ±0.07. The large- 
scale quasar linear bias is estimated to be bQ = 1.38 ± 
0.10 at (z) ~ 0.53. This bias infers that quasars at 
these redshift reside in halos with typical mass of ~ 
4 x 10 12 /j _1 M q (using the Tinker et al. 2005 fitting 
formula), similar to quasar clustering measurements 
at high-redshift, but lower than the typical halo mass 
~ 10 13 Ii~ 1 Mq for massive galaxies in SDSS. Thus most 
of these low-redshift quasars are not the descendants 
of their high-redshift counterparts, which would have 
evolved into more massive and more biased systems 
(such as the hosts of CMASS galaxies). 

• We found weak luminosity dependence of the large- 
scale quasar linear bias, over the luminosity range 
-23.5 > Mi(z = 2) > -25.5 probed by our sample. 
This result is generally consistent with other quasar 
clustering measurements at different redshifts. This 
weak luminosity dependence suggests that quasars with 
fixed luminosity spread over a broad range of host 
halo masses, in qualitative and quantitative agreement 
with predictions from several theoretical models (e .g., 
ILidz et dJ200alSr^l200^IConrov & Whitdl2013l) . 

• We performed HOD and mock catalog-based model- 
ing of the measured CCF. For the HOD modeling, we 
found large degeneracies in the HOD parameteriza- 
tions such that different HODs can reproduce the CCF 
equally well, with different host halo mass distributions 
and satellite fractions. This result highlights the limita- 
tions and ambiguities in the standard HOD approach 
for modeling the quasar population. Additional infor- 
mation is needed in order to break the degeneracies in 
the quasar HOD models. 

For the mock-based approach, we found the simple 
model in lConrov & White! d2013l) that relates quasars to 



galaxies can reproduce the CCF reasonably well. Un- 
der such a model framework, we need a satellite frac- 
tion of quasars (i.e., fraction of quasars hosted by satel- 
lite galaxies) of / sat ~ 10%. Just as for the HOD-based 
modeling, however, we cannot rule out other models by 
which quasars can inhabit dark matter halos and pro- 
duce the same CCF. 

The difficulty of finding a unique HOD model for 
quasars probably lies primarily in the fact that quasars 
are a sparse population with an unknown duty cycle rel- 
ative to halos (or galaxies). The large scatter between 
quasar luminosity and halo mass also makes it difficult 
to use luminosity-dependent clustering as an additional 
constraint in quasar HOD modeling. 

With the upcoming data release of the BOSS survey, we 
will eventually have a spectroscopic CCF sample with ~ 50% 
more quasars and more CMASS galaxies with the final SDSS- 
III data release. The new data will increase the cross-pair 
counts by ~ 50%. On small scales (r p < l/i _1 Mpc) where 
Poisson statistics dominate, we therefore expect ~ 20% im- 
provement in the errors of w p measurements. These changes 
will potentially be able to reveal differences in the small-scale 
clustering when binned in quasar luminosity. In the short 
term, we also plan to measure the CCF using spectroscopic 
SDSS-DR7 quasars and the photometric CMASS galaxy sam- 
ple, which will have the same cross-sample coverage as the 
final spectroscopic CMASS sample and is free of fiber col- 
lision losses. Future deeper galaxy and quasar surveys over 
large areas can improve the pair statistics further, and at the 
same time increase the dynamical range in quasar luminosity. 
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Figure 1. Correlation matrix of w p (r p ) for the full sample cross correlation 
(DR10 CMASS galaxies with DR7 uniform quasars). This is the normalized 
covariance matrix, i.e. correlation matrix, such that the diagonal elements are 
unity calculated on 50 jackknife samples. The contours correspond to values 
of 0.75, 0.50, and 0.25. 
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APPENDIX 

We estimate errors on our clustering measurements us- 
ing the jackknife resampling technique (as discussed in Sec- 
tion [3}. We use the full covariance, which includes the cor- 
relation between bins in the correlation function as shown in 
Figure Q] We use a fiducial value of 50 jackknife regions, 
which we define such that each region has the same unmasked 
area on the sky and is roughly rectangular (where possible). In 
this section, we evaluate some of the effects on the errors due 
to varying the number of jackknife regions for our measure- 
ments of the projected two-point correlation function. Specif- 
ically, we compare error estimates on our cross correlation 



measurement for the full sample using 10, 25, 50, 75, and 100 
jackknife regions. 

The number of jackk nife regions is somew hat arbitrary (see 
detailed discussion in Norberg et al. 2009). Using too few 
jackknife samples will result in a low number of realizations 
to estimate the variance, and can formally cause the covari- 
ance matrix to become singular (when the number of samples 
is less than the number of bins). The use of too many jack- 
knife regions causes each region to become small in area (and 
therefore volume) and can inaccurately represent the cosmic 
(sample) variance in the large-scale errors. At a minimum, we 
must ensure the size of each jackknife region is significantly 
larger than the largest scales we measure in the data. 

We first investigate the magnitude of the diagonal errors, 
which we show in Figure|2] The values of a can vary by up to 
30-40%, but are otherwise roughly equivalent. There is no 
systematic bias in the values that affects one choice more than 
any other across all the bins. A lower number of jackknife 
samples, however, results in a larger variation in the values, 
as we would expect. 

To quantify how well we resolve the structure of the cor- 
relation matrix (e.g. Figure Q3, we perform a singular value 
decomposition (SVD) on the correlation matrix. The SVD 
effectively rotates the matrix into an orthogonal space which 
can be thought of as a combination of eigenvectors and eigen- 
values. The singular values (SVs) are eigenvalues (defined 
to be positive) which are the multiplicative amplitude of the 
corresponding (normalized) eigenvector. The SVs are typi- 
cally numbered such that they are monotonically decreasing, 
and can be interpreted as a measure of the "importance" of 
each mode in terms of contributing to the observed structure 
in the full correlation matrix. For example, an N by N diag- 
onal correlation matrix (i.e. the identity matrix) would result 
in N SVs that were all equal in value. The ratio of the largest 
SV divided by the smallest SV is referred to as the condition 
number, and if significantly large can result in poor numeri- 
cal results when the matrix is inverted (i.e. an ill-conditioned 
matrix) which is performed in model fitting. 

We show the SVs for our correlation matrices in Figure 
We clearly see our expectation of the ill-conditioned matrix 
for Afj ac k =10 since we are using 22 bins. A larger N-^y re- 
sults in a better conditioned matrix (a line that appears more 
flat as the SVs vary less). We also notice quickly diminishing 
returns for larger numbers of samples: while there is a dra- 
matic difference between 10 and 50 samples, it is much less 
of a difference for the larger numbers of jacknife samples. 

We conclude from these investigations that using less than 
50 jackknife samples could be troubling. Taking into account 
the area coverage of our data (4122 deg 2 ), 50 jackknife sam- 
ples result in each jackknife region covering about 82 deg 2 
(roughly 9 deg or less on a side). As our statistical errors are 
signifi cantly larger than the galaxy autocorrelation function 
(e.g. IZehavi et alj|201 ll) . we are not overly concerned about 
resolving each element of the covariance matrix. Our choice 
of 50 jackknife samples is a factor of 2 larger than the number 
of bins. 
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Table 1 

Normalized covariance matrices from the full CCF sample and from individual CCF subsamples. The corresponding diagonal errors are tabulated in Table|2] A 
portion (i.e., the covariance matrix from the full CCF sample) is shown here for its content. The table is available in its entirety in the electronic version of this 



paper. 


r p |/r'Mpc| 


0415 


0.154 


0.205 


0.274 


0.365 


0.487 


0.649 


0.866 


1.155 


1.540 


2.054 


2.738 


3.652 


4.870 


6.494 


8.660 


1 1 .548 


15.399 


20.535 


27.384 


36.517 


48.697 


0.1 15 


1 .000 


-0.002 


-0.087 


-0.047 


0.001 


0.137 


0.015 


0.141 


-0.019 


-0.023 


-0.161 


-0.035 


-0.083 


0.273 


0.305 


0.049 


0.114 


-0.086 


0.078 


0.123 


0.022 


-0.182 


0.154 




1.000 


0.082 


0.106 


0.028 


-0.010 


-0.179 


0.016 


0.058 


-0.029 


-0.116 


0.104 


-0.215 


0.000 


-0.077 


0.197 


0.008 


0.147 


0.125 


0.147 


0.107 


0.136 


0.205 






1 .000 


0.055 


-0.097 


-0.042 


0.230 


0.126 


-0.164 


-0.127 


-0.001 


0.101 


0.057 


-0.155 


-0.220 


0.013 


-0.015 


0.109 


-0.070 


-0.069 


0.091 


0.052 


0.274 








! .000 


0.226 


0.058 


-0.375 


-0.055 


0.074 


-0.213 


0.021 


-0.078 


-0.050 


0.034 


-0.061 


0.229 


0.229 


0.389 


0.262 


0.248 


0.213 


0.087 


0.365 










1.000 


0.071 


-0.118 


0.340 


0.220 


-0.050 


0.138 


0.284 


0.034 


0.110 


-0.006 


0.06S 


0.089 


0.144 


0.002 


-0.067 


-0.259 


0.030 


0.487 












1.000 


0.027 


-0.049 


-0.061 


0.09 1 


-0.068 


-0.172 


-0.072 


-0.077 


-0.253 


0.047 


-0.038 


0.041 


-0.190 


-0.190 


-0.214 


-0.302 


0.649 














1.000 


0.154 


-0.070 


0.010 


0.309 


0.133 


0.019 


-0.007 


-0.087 


0.098 


0.053 


-0.094 


-0.068 


-0.172 


-0.189 


-0.210 


0.S66 
















1.000 


0.024 


0.033 


0.266 


0.269 


-0.108 


0.144 


0.150 


0.147 


0.051 


0.092 


0.029 


0.158 


-0.088 


-0.068 


1.155 


















1.000 


-0.087 


0.140 


0.180 


0.070 


0.170 


0.389 


0.123 


0.124 


0.194 


0.082 


-0.064 


-0.022 


0.133 


1.540 




















1 .000 


-0.087 


0.031 


0.105 


0.282 


0.276 


-0.060 


0.084 


0.020 


0.190 


0.333 


0.233 


0.157 


2.054 






















1.000 


0.231 


0.047 


0.091 


0.097 


0.076 


0.2S0 


0.245 


0.127 


0.072 


-0.137 


-0.165 


2.738 
























1.000 


0.389 


0.350 


0.166 


0.297 


0.153 


0.336 


0.285 


0.077 


-0.149 


0.1 1 I 


3.652 


























1 .000 


0.110 


0.269 


0.173 


0.175 


0.216 


0.262 


0.063 


-0.011 


0.208 


4.870 




























1 .000 


0.487 


0.429 


0.268 


0.245 


0.355 


0.270 


0.1 19 


0.004 


6.494 






























1 .000 


0.363 


0.346 


0.311 


0.407 


0.423 


0.275 


0.236 


8.660 
































1.000 


0.496 


0.598 


0.547 


0.294 


0.224 


0.207 


11.548 


































1.000 


(1.601 


0.565 


0.425 


0.327 


(1.1X1 


15.399 




































1 .000 


0.703 


0.548 


0.440 


0.385 


20.535 






































1.000 


0.679 


0.488 


0.482 


27.384 








































1 .000 


0.649 


0.403 


36.517 










































1.000 


0.694 


48.697 












































1.000 



Gunn, J. E., et al. 1998, AJ, 1 16, 3040 
— . 2006, AJ, 131,2332 

Guo, H., Zehavi, I., & Zheng, Z. 2012a, ApJ, 756, 127 
Guo, H., et al. 2012b, ApJ, submitted, arXiv:1212.1211 
Haiman, Z., & Hui, L. 2001, ApJ, 547, 27 
Hennawi, J. R, et al. 2006, AJ, 131, 1 
Hickox, R. C, et al. 2009, ApJ, 696, 891 
—.2011, ApJ, 731, 117 

Hirschmann M., Somerville R. S., Naab T, Burkert A. 2012, MNRAS, 426, 
237 

Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 2001, AJ, 122, 
2129 

Hopkins, R F, Hernquist, L., Cox, T. J., & Keres, D. 2008, ApJS, 175, 356 
Hopkins, R F, Hernquist, L., Martini, R, Cox, T. J., Robertson, B., Di 

Matteo, T, & Springel, V. 2005, ApJ, 625, L71 
Hopkins, R F, Lidz, A., Hernquist, L., Coil, A. L., Myers, A. D., Cox, T. J., 

& Spergel, D. N. 2007a, ApJ, 662, 110 
Ivashchenko, G, Zhdanov, V. I., & Tugay, A. V. 2010, MNRAS, 409, 1691 
Ivezic, Z, et al. 2004, Astronomische Nachrichten, 325, 583 
Kaiser, N. 1987, MNRAS, 227, 1 
Kayo, I. and Oguri, M. 2012, MNRAS, 424, 1363 
Komatsu, E., et al. 2011, ApJS, 192, 18 
Krumpe, M., Miyaji, T, & Coil, A. L. 2010, ApJ, 713, 558 
Krumpe, M., Miyaji, T, Coil, A. L., & Aceves, H. 2012, ApJ, 746, 1 
Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 

Li, C, Kauffmann, G., Wang, L., White, S. D. M., Heckman, T. M., & Jing, 

Y. P. 2006, MNRAS, 373, 457 
Lidz, A., Hopkins, P. R, Cox, T. J., Hernquist, L., & Robertson, B. 2006, 

ApJ, 641,41 

Lupton, R., Gunn, J. E., Ivezic, Z., Knapp, G. R., & Kent, S. 2001, in 
Astronomical Society of the Pacific Conference Series, Vol. 238, 
Astronomical Data Analysis Software and Systems X, ed. F. R. Harnden, 
Jr., F. A. Primini, & H. E. Payne, 269 

Maraston, C, et al. 2012, arXiv: 1207.61 14 

Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Martini, P., & Weinberg, D. H. 2001, ApJ, 547, 12 

McBride, C. K. and Connolly, A. J. and Gardner, J. P. and Scranton, R. and 

Newman, J. A. and Scoccimarro, R. and Zehavi, I. and Schneider, D. P. 

2011, ApJ, 726, 13 
Miyaji, T, Krumpe, M., Coil, A. L., & Aceves, H. 201 1, ApJ, 726, 83 
Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 
Myers, A. D., Branner, R. J., Nichol, R. C, Richards, G. T, Schneider, 

D. P., & Bahcall, N. A. 2007a, ApJ, 658, 85 
Myers, A. D., Branner, R. J., Richards, G. T, Nichol, R. C, Schneider, 

D. P., & Bahcall, N. A. 2007b, ApJ, 658, 99 
Myers, A. D., Richards, G. T, Branner, R. J., Schneider, D. P., Strand, N. E., 

Hall, P. B., Blomquist, J. A., & York, D. G. 2008, ApJ, 678, 635 
Myers, A. D., et al. 2006, ApJ, 638, 622 

Navarro, J. R, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Norberg, P., et al., 2001, MNRAS, 328, 64 

Norberg, P., Baugh, C. M., Gaztafiaga, E., & Croton, D. J. 2009, MNRAS, 

396, 19 

Nuza, S. E., et al. 2012, MNRAS, in press, arXiv: 1202.6057 
Padmanabhan, N., White, M., Norberg, P., & Porciani, C. 2009, MNRAS, 
397, 1862 

Parejko, J., et al. 2012, MNRAS, submitted 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, S. M., 

Lupton, R. H, & Ivezic, Z. 2003, AJ, 125, 1559 
Porciani, C, Magliocchetti, M., & Norberg, P. 2004, MNRAS, 355, 1010 
Porciani, C, & Norberg, P. 2006, MNRAS, 371, 1824 



Richards, G. T, et al. 2002a, AJ, 123, 2945 
Richards, G. T. et al. 2006, AJ, 131, 2766 

Richardson, J., Zheng, Z., Chatterjee, S., Nagai, D., & Shen, Y. 2012, ApJ, 
755, 30 

Ross, N. P., et al. 2009, ApJ, 697, 1634 

Ross, N. P., et al. 2012, ApJS, 199, 3 

Sanchez A.G., et al., 2012, submitted, arXiv: 1203.6616 

Schlegel D., White M., Eisenstein D.J., 2009, The Astronomy and 

Astrophysics Decadal Survey, Science White Papers #314, 

arXiv:0902.4680 
Scranton, R., et al. 2002, ApJ, 579, 48 
Schneider, D. P., et al. 2010, AJ, 139, 2360 

Shankar, F, Weinberg, D. H, & Miralda-Escude, J. 2009, ApJ, 690, 20 
Shankar, F, Crocce, M., Miralda-Escude, J., Fosalba, P., & Weinberg, D. H. 

2010, ApJ, 718, 231 
Shankar, F, Weinberg, D. H, & Shen, Y. 2010, MNRAS, 406, 1959 
Shanks, T, Croom, S. M., Fine, S., Ross, N. P., & Sawangwit, U. 201 1, 

MNRAS, 416, 650 
Shaver, P. A. 1984, A&A, 136, L9 
Shen, Y. 2009, ApJ, 704, 89 

Shen, Y, Greene, J. E., Strauss, M. A., Richards, G. T, & Schneider, D. P. 

2008, ApJ, 680, 169 

Shen, Y, Strauss, M. A., Hall, P. B., Schneider, D. P., York, D. G., & 

Bahcall, N. A. 2008, ApJ, 677, 858 
Shen, Y, & Kelly, B. C. 2010, ApJ, 713, 41 
— . 2012, ApJ, 746, 169 
Shen, Y, et al. 2007, AJ, 133, 2222 
— . 2009, ApJ, 697, 1656 
— . 2010, ApJ, 719, 1693 
—.2011, ApJS, 194,45 

Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 

Silk, J., & Rees, M. J. 1998, A&A, 331.L1 

Smee, S. A., et al. 2012, AJ, submitted, arXiv: 1208.2233 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Tegmark, M., & Peebles, P. J. E. 1998, ApJ, 500, L79 

Thacker, R. J., Scannapieco, E., Couchman, H. M. P., & Richardson, M. 

2009, ApJ, 693, 552 

Tinker, J. L., Weinberg, D. H, Zheng, Z., & Zehavi, I. 2005, ApJ, 631, 41 
Tucker, D. L., et al. 2006, Astronomische Nachrichten, 327, 821 
van den Bosch, R, More, S., Cacciato, M., Mo, H, & Yang, X. 2012, 

arXiv: 1206.6890 
Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 
Wake, D. A., Croom, S. M., Sadler, E. M., & Johnston, H. M. 2008, 

MNRAS, 391, 1674 
White M., 2002, ApJS, 579, 16 

White, M., Martini, P., & Cohn, J. D. 2008, MNRAS, 390, 1 179 
White, M., Zheng, Z., Brown, M.J.I., Dey, A., Jannuzi, B.T., 2007, ApJ, 655, 
L69 

White, M., et al. 201 1, ApJ, 728, 126 

White, M., et al. 2012, MNRAS, 424, 933 

Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 

York, D. G., et al. 2000, AJ, 120, 1579 

Yu, Q., & Lu, Y. 2004, ApJ, 602, 603 

— . 2008, ApJ, 689, 732 

Zehavi, I., et al. 2005, ApJ, 630, 1 

—.2011, ApJ, 736,59 

Zheng, Z. 2004, ApJ, 610, 61 

Zheng, Z., Berlind, A. A., Weinberg, D. H, et al. 2005, ApJ, 633, 791 
Zheng, Z., Coil, A. L., & Zehavi, I. 2007, ApJ, 667, 760 



